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Abstract 

In this paper we consider a large set of variable coefficient linear 
systems of ordinary differential equations which possess two different time 
scales, a slow one and a fast one. A small parameter e characterizes the 
stiffness of these systems. We approximate a system of o.d.e.s in this set by 
a general class of multistep discretizations which includes both one-leg and 
linear multistep methods. We determine sufficient conditions under which each 
solution of a multistep method is uniformly bounded, with a bound which is 
independent of the stiffness of the system of o.d.e.s, when the step size 
resolves the slow time scale but not the fast one. We call this property 
stability with large step sizes. 

The theory presented in this paper lets us compare properties of one-leg 
methods and linear multistep methods when they approximate variable 
coefficient systems of stiff o.d.e.s. In particular, we show that one-leg 
methods have better stability properties with large step sizes than their 
linear multistep counterparts. This observation is consistent with results 
obtained by Dahlquist and Lindberg [11], Nevanlinna and Liniger [32] and van 
Veldhuizen [41]. Our theory also allows us to relate the concept of D- 
stability (van Veldhuizen [41]) to the usual notions of stability and 
stability domains and to the propagation of errors for multistep methods which 
use large step sizes. 
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1 . INTRODUCTION 


Model equations (usually obtained through a linearization procedure) have 
played an extremely important role in assessing the behavior of a 
discretization of a system of stiff differential equations. The test equation 

(1.1) = * e $ is a constant, 

has led to the concepts of A-stabllity [7], A(ct)-stability [43], stiff- 
stability [17], and L-stability [13], However, this model is too simple. 
Examples due to Gourlay [18] show that two different one-step methods with the 
same domain of absolute stability can have entirely different properties when 
they approximate the variable coefficient test equation 

(1.2) g.= X (t)y, X(t)€$ for 0 < t < T. 

For this reason and others Prothero and Robinson [35] introduced a theory 
for the stability and accuracy of one-step discretizations of stiff o.d.e.s 
based on the model equation 

(1.3) jfc = g'(t) + X(y - g(t)) 

where X is a complex constant with negative real part, ( )' denotes the 
derivative with respect to t, and g'(t) is any function that is defined and 
bounded for t € [0,T] where T is some constant. 

In several papers nonlinear systems of o.d.e.s which satisfy a 
monotonicity condition have been used for model problems. See Dahlquist [8], 
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Odeh and Liniger [34], Stetter [37], Burrage [3], Burrage and Butcher [4], 
Butcher [5], and Crouzeix [6], 

Most stiff systems, for example the equations which arise in chemistry or 
electric circuit theory, have several components which respond with widely 
differing time constants. Furthermore, these components are usually coupled, 
and both the coupling and the time constants are time dependent. (cf. Bjurel, 
et al., [2], Lapidus and Schiesser [24] or Willoughby [44].) For these 
reasons models (1.1) - (1.3) and equations with monotone nonlinearities are 
still inadequate, and several authors including Dahlquist [9], Miranker [30], 
[31], van Veldhuizen [39], [40], [41], and Stetter [36] have proposed more 

refined model equations for discretizations of stiff o.d.e.s. Stetter 
suggested that an appropriate model equation should possess the following 
properties : 

1. It should permit the simultaneous occurrence of slowly varying and 
rapidly varying solution components. 

2. It should have a Jacobian matrix with a time-dependent eigensystem. 

3. It should contain a small parameter to permit the consideration of a 
limit process corresponding to a transition to arbitrarily high 
stiffness . 

In this paper we will consider model linear systems of o.d.e.s which 
satisfy properties 1-3. These equations are given by 

(1.4) = A(t,e)v, 0 < e £ e 0 , v(0) given, 

and 

= A(t,e)y + F(t,e), 0 < t < T, 0 < e £ e Q , y(0) given, 


(1.5) 
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where A(t,e) is a n x n matrix, v, y and F(t,e) are n-component 
vectors, and T, are positive constants with « 1. A(t,e) and 

F(t,e) should satisfy the conditions in the following three assumptions: 


Assumption 1.1. There exists an invertible matrix T(t,e) with T(t,e), 
T ^(t,e)£C P (t,e,B)^) and p _> 3 so that 


( 1 . 6 ) 


T 1 (t ,e )A( t ,e )T( t ,e ) 


D, i (t,e) 

— 0 

e 

0 D 22 (t,e) 


D ( t , £ ) • 


Here D^(t,e) € C p (t ,e ,B) is an m x m matrix with 1 _< m < n and 
D 22 (t,e) £ C^(t,e ,B) is an (n-m) x (n-m) matrix. 

Assumption 1.2. D^(t,e) is an invertible matrix for 0 _< t _< T and 

0 <_ e <_ eg with j(t,e) £ C p (t,e,B) . Furthermore, the eigenvalues of 
D^Ct.e), denoted by X i (t,e), i = l,««»,m, satisfy Re{x i (t,e)} < 0, 

1 = 1 , • * * ,m, and are distinct for 0 < t _< T and 0 _< e _< gg . 

Assumption 1.3. T (t,e)F(t,e) = (— f ^(t,e) ,f ^.j.(t,e)) where 
f-j.(t,e) £ C p (t,e,B) denotes an m -component vector and f^Ctje) £ C P (t,e,B) 
denotes an ( n-m ) -component vector. 


( 1 ) 


All notation is defined in Part I of the Appendix. 
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An asymptotic decomposition for the solutions of systems (1.5) is easily 
obtained. The results in Hoppensteadt [19], Hoppensteadt and Miranker [20], 
and Kreiss [22] imply that if e is sufficiently small, then the solutions of 
systems (1.5) which satisfy the conditions in Assumptions 1.1 - 1.3 have an 
additive decomposition with the form 


(1.7) 


s R 

y(t ,e) = y (t,e) + y (x,e) with t = t/e. 


y (t,e) is a slowly varying solution of (1.5) with (p - 2) derivatives with 

R 

respect to t bounded independently of e. y (x,e) is a rapidly varying 
solution of (1.5) which is highly oscillatory (rapidly decaying) if 
Re{X i (t,e)} = 0, i = l,***,m, (Rejx^t.e)} < 0, i = l,»**,m). Both y S (t,e) 
and y^(x,e) are uniquely determined up to terms 0(e^ 

In order to compare the properties of various multistep methods, we will 
consider the general class of multistep approximations of (1.4) and (1.5) 
given by 


(1.8) £ (a. I + kb.A(x.(t ),e))v n 1 = 0, n = r,»*»,N, 

i=0 1 


(i.9) l (a i I + kb i A(T i (t n ),e))y n + k F (T Q (t n ) , • • • ,x r (t n ) ,k,e) = 0, 
i=0 

n = r , • • • ,N. 


Here a^, b-p i = 0,»»»,r, denote given constants, k > 0 denotes a constant 
step size, t n = nk, n = 0,1, x^t) € C^(t,B) , i = 0,»**,r, x^t) = 

x Q (t) + 0(k) for i = !,•••, r, and 0 < t < T, 0 <_ x i (t) < T for 
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i = 1, and t^ _< t £ T, and F(-Tg(t) , • • • ,T r (t) ,k,e) denotes a function 
which depends on F(x^(t),e), i = 0,»»»,r. We will assume that F = 0 when 
F(t,e) = 0 and that the coefficients in (1.8) and (1.9) are normalized so 
that ag = 1. We will refer to method (1.8) as the homogeneous multistep 
method corresponding to method (1.9). 

Under special conditions methods (1.9) lead to well-known classes of 

raultistep methods. If we set x^(t) = t-ik, i = 0,»»»,r, and 
r 

F= £ b.F(t-ik,e), then (1.9) becomes 
i=0 1 

r r 

(1.10) Y (a. I + kb.A(t .,e))y n 1 + k T b.F(t .,e) = 0 

. L . K i l n-i’ . L n i n-i’ 

i=0 i=0 

which is a linear multistep method. On the other hand, if we set b^ = 1 , 

r 

i = (),•••, r, t . ( t ) = T n (t) = y S.(t-jk), i = 0,«**,r, where the g.'s are 
1 j=0 J J 

given constants and F= F(Tg(t),e), then (1.9) becomes 

(1.11) T (a. I + kA/ y e.t . .eUy 11 " 1 + kF / T g.t .,e\ = 0 

i=0 \ 1 \j=0 J // \j=0 n_ '* / 

which is a one-leg method. These methods were introduced by Dahlquist in [8]. 

In this paper we are particularly interested in the properties of 
discretizations (1.8) and (1.9) when the step size k resolves only the slow 
time scales in system (1.4) or (1.5) and not the fast ones. This is the 
typical situation for discretizations of stiff equations. Therefore, we are 
interested in the cases when e tends to zero for fixed k > 0 or k 0 
with k/e bounded away from zero. (The limit process k/e 0 would take us 

back to the well-known Dahlquist-Henrici theory of discretizations on 
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"suf f Iciently fine" grids.) Consequently, the two parameters e and k will 
vary in the trapezoidal region 


(1.12) r = {(k,e)€]R?:0<e^eQ and a e _< k _< k Q } 

k 

where a and kg are two positive constants (see Fig. 1). 



Following the ideas and results in van Veldhuizen [39], [40], Stetter [36] 
suggested that a desirable discretization of a system (1.5) should satisfy two 
properties. First, each solution of the numerical method should be uniformly 
bounded with a bound which is independent of the stiffness of the equation. 
We formulate this requirement as 

Property 1.1 . For some trapezoidal region r, each solution of (1.9), 
denoted by y n , n = Q,.«*,N, should be uniformly bounded for all (k,e) € T. 
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A multistep approximation of a system (1.5) which satisfies Property 1.1 

is called a stable approximation of (1.5) with large step sizes. 

The second property pertains to the accuracy of the numerical solution. 

In particular, only the slow time scales in (1.5) are resolved by the step 

size while the fast scales are not resolved. Therefore, the numerical 

approximation of (1.5) should have a two-term decomposition corresponding to 

(1.7). The first terra should approximate the slowly varying solution 

y S (t,e) in the decomposition (1.7) with some degree of accuracy. The second 

term should decay with time. This behavior of the second terra is consistent 

R R 

with the properties of y (t,e) when y (t,e) is rapidly decaying. In the 
case when y (x,e) contains rapidly oscillating components, this restriction 
assumes that the rapidly oscillating components have no interest and should be 
attenuated. We formulate this accuracy requirement as 

Property 1.2 . For each value of n = 0,1,»«*,N, y n should have the 
decompositon 

(1.13) y 11 = y 11 + E n 

s 

n s 

where y^ approximates y (t^,£) uniformly well for all (k,e) £ T, that is, 

(1.14) lly^ - y S (t n ,e)|| < (k + e ) when (k,e)£ T 

where Kj is some constant independent of k and e and q^, q2 are two 

positive integers. E n should decay, that is, 

II E n II _< K^S 11 when (k,c) €T 


(1.15) 
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where K.£ and 0 < 5 < 1 are constants independent of k and e. 

In this paper we apply the concept of stability domains together with the 
usual Dahlquist-Henrici stability theory to determine sufficient conditions 
under which a multistep approximation of (1.5) with the form (1.9) satisfies 
Property 1.1. One immediate consequence of this theory is that one-leg 
methods have better stability properties with large step sizes than their 
linear multistep counterparts when they approximate variable coefficient (and 
nonlinear) systems of stiff ordinary differential equations. A second 
consequence of our theory is to relate the concept of D-stability (van 
Veldhuizen [41]) to the usual notions of stability and stability domains and 
to the propagation of errors for multistep methods which use large step 
sizes. Additional consequences are discussed in Sections 3 and 5. We 
determine sufficient conditions under which a multistep approximation of (1.5) 
satisfies Property 1.2 in a forthcoming paper using the techniques developed 
in [28]. 

The results in this paper either extend or complement known results in 
several ways. Van Veldhuizen [39], [40], established Properties 1.1 and 1.2 
for a class of one-step methods using techniques which are entirely different 
from the ones presented in this paper. He was unable to extend his theory to 
multistep methods. Dahlquist and Lindberg [11] compared the stability, 
accuracy, and step changing properties of the implicit midpoint rule and the 
trapezoidal rule when these methods approximate stiff linear and nonlinear 
systems of o.d.e.s. Kreiss [21] essentially determines sufficient conditions 
under which A(a)-stable linear multistep approximations of systems (1.4) 
satisfy Property 1.1 when the matrix A(t,e) in (1.4) is of essentially 
diagonally dominant type. Kreiss's results cannot be extended to the class of 



-9- 


systems (1.4) which satisfy Assumptions 1.1 and 1.2. (See van Veldhuizen [41] 
for details.) Dahlquist and Soderland [12] essentially establish Property 1.1 
for a class of multistep and multistep compound discretizations of systems 
(1.5) which appear in singular perturbation form using the theory of G- 
contractivity. The systems (1.5) which we consider do not appear in singular 


perturbation form. 

Furthermore, 

we consider 

a class 

of 

multistep 

discretizations of 

(1.5) which is 

larger than 

the class 

of 

multistep 


discretizations considered by Dahlquist and Soderland. Finally, Abrahamson, 
Keller, and Kreiss [1] have established Properties 1.1 and 1.2 for a class of 
three-point difference approximations of some linear two-point boundary value 
problems which appear in singular perturbation form. 

This paper is organized in the following way: In Section 2 we consider 

multistep approximations of the special systems s = D^(t,e)y where 
Dn(t,e) satisfies the conditions in Assumption 1.2. These systems govern 
the rapidly varying components of systems (1.4) and (1.5). We use properties 
of stability domains to determine a trapezoidal region r with the form 
(1.12) so that all solutions of a multistep approximation of these special 
systems decay with time when (k,e) £T. In Section 3 we determine sufficient 
conditions under which a multistep approximation of system (1.5) satisfies 
Property 1.1. We also present a thorough discussion of the practical 
implications of Property 1.1 for a multistep method. In Section 4 we relate 
the concept of D-stability (van Veldhuizen [41]) to the results obtained in 
Section 3. In Section 5 we present some numerical results. 
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2. MULTISTEP APPROXIMATIONS OF THE UNRESOLVED FAST SCALES 

In this section we will consider multistep approximations of the special 
systems 

(2.1) e = D 11 (t,e)w, 0_<t_<T, 0 < e £ e Q , w(0) given, 

where T and « 1 are given constants and D^(t,e) £C^(t,e,B) with 
p > 1 is an m x m matrix which satisfies the conditions in Assumption 
1.2. Systems (2.1) govern the rapidly varying components of systems (1.4) and 
(1.5). In particular the results in [16] or [42] show that if £q is 
sufficiently small, then for any integer 0 < v £ p-1 the solutions of (2.1) 
have an asymptotic decomposition with the form 

/ / S 1 (y,s)dY / 8 m (Y»e)dY\ 

w(t) = diag \e^ ,»»»,e^ / w(0) + 0(e V ) 

for 0 < t < T and 0 < e £ £q 


where 


X.(t,0) 

3.(t,e) = 


,m, 0 £ t £ T and 0 < e £ Sq. 


e 


+ 0(1) for i = 1 , • • • 
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Consequently, the solutions of (2.1) are highly oscillatory (rapidly decaying) 
if Re{X i (t,e)} = 0, i = l,***m, (Re{X i (t,e)} < 0, i = l,«»»,m). 

We will consider multistep approximations of (2.1) with the form (1.8) 
using a fixed step size k which is bounded but does not resolve the rapidly 
varying components, that is, k and e satisfy k _< k^ and k/e » 1 
where kg is some constant. We use properties of stability domains to 
determine sufficient conditions under which all solutions of a multistep 
approximation of (2.1) decay with time when k and e satisfy the previously 
stated conditions. This behavior of the solution of the multistep method is 
consistent with the properties of the solutions of (2.1) when the solutions of 
(2.1) decay. In the case when the solutions of (2.1) contain rapidly 
oscillating components, this restriction assumes that the rapidly oscillating 
components have no interest and should be attenuated. 

In order to state our main result, we need to define the usual notion of a 
domain of absolute stability for a multistep method (1.9) and to introduce two 
concepts which generalize the notions of A-stability and L-stability to 
subsets of the complex plane other than the left half plane. 

Approximate the test equation (1.1) by a multistep method (1.9) and set 
z = kX to obtain the difference equation 


( 2 . 2 ) 


r 

1 (a, + b z)w n 1 = 0, r _< n < °°. 

i=0 1 1 


Set 


~n , n 
w = (w , 


n-r+1 J 
>w ) 


for r-1 < n < ® and let ft denote a subset of 


the complex plane 
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Def Inition 2.1. 

a) The domain of absolute stability of (1.9) consists of all 

points such that (ag + bgz) - ^ exists and there exists a constant 

K(z) so that all solutions of (2.2) satisfy 

l!w n || _< K(z)llw~*l| for r-1 _< j <_ n <_ °°. 

b) The multistep method (1.9) is strictly fl-stable if (ag + bgz) - ^ 
exists for all zGfl, and for each z e 51 there exist constants K(z) and 
6(z) with 0 _< 6(z) <1 so that all solutions of (2.2) satisfy 

11 w n It < K(z)(5(z)) n ^ II II for r-1 < j £ n < ». 

c) The multistep method (1.9) is strongly ^-stable if 

(2.3) (ag + bgZ ) * exists and is uniformly bounded for all z € Cl, 

and there exist constants K > 0 and 0 _< 6 < 1, independent of z, so that 
all solutions of (2.2) satisfy 

llw n || _< K6^ n ^Hw^l for r-1 _< j _< n < « for all z € Cl, 

Algebraic conditions which determine the domain of absolute stability of a 
multistep method or subsets Cl in which a multistep method is either strictly 
or strongly n-stable are easily stated. 
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The solutions of the difference equation (2.2) are determined by the roots 

r _ . 

of the polynomial R(z,k) = £ (a. + b.z)< which we denote by k.(z), 

i=0 1 1 1 ' 

i = 1 , • • • ,r. 

Definition 2.2 . The roots of a polynomial q(x) satisfy the root 
condition if all roots of q(x) = 0, denoted by x^, v = 1, •••,(), satisfy 
l x v l 1 » v ** 1 , • • • ,Q, and those roots which lie on the unit circle are 
simple. 

Lemma 2.1. 

a) The domain of absolute stability for method (1.9) consists of all 

points z £ (£ such that (ag + bgz) ^ exists and the roots of R(z,k) 

satisfy the root condition. 

b) The multistep method (1.9) is strictly ft -stable if (ag + bgz) ^ 

exists for all z € ft and the roots of R(z,k) satisfy |tc^(z)| < 1, 

i = 1 , • • • ,r. 

c) The multistep method (1.9) is strongly ft -stable if condition (2.3) is 
satisfied and there exists a constant 0 _< y < 1, independent of z, so that 
the roots of R(z,k) satisfy | ( z ) | < 6 for i = l,...r, and all z € ft. 

Example 2.1. Consider the trapezoidal approximation of (1.5) given by 

y n = y n 1 + j ( A (t n ,e)y n + A (t n _ 1 ,e)y n_1 ) + (F(t n , e ) + F(t n _i> e) )* 


Set D = { (x + iy) € <£ : x _< 0} and D' = { (x + iy) € £ : x < 0} . Then 
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i) The domain of absolute stability of the trapezoidal rule is the set D. 
ii) The trapezoidal rule is strictly D'-stable but not strongly D'-stable. 
Identical remarks apply to the implicit midpoint rule 


n n-1 , k n , n-K , , „ 

y =y + J A(t n - V? ’ e)(y 7 > + kF(t n - V? ,e) 


which is the one-leg twin associated with the trapezoidal rule. 

Example 2.2 . For each value of p = the backward difference 

formula of order p ([23], p. 242 or [17], Chapter 11) is strongly 
ftp-stable in the region 


(2.4) ftp = {z e<t : -a + e _< tt - arg(z) _< a p - e and |z| > p} 


where p is any positive constant and e is any constant satisfying 
0 < e <_ a with = 90°, 02 = 90°, = 88°, = 73°, = 51° and 


If we set p = 0 in (2.4), then each backward difference formula of 
order p is strictly ftp-stable but not strongly ftp-stable. 

We are now ready to state the main result of this section. 


Theorem 2.1. Approximate system (2.1) by a multistep method (1.8) with a 


step size k satisfying 0 < k < k where k is some constant and set 


~n , n 

y = (y , 


H” X T 

•••,y ) . Let ft denote a subset of the complex plane. Assume 


that the multistep method is implicit, strongly ft-stable^ and there exists a 


trapezoidal region Y' with the form 
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(2.5) r' = {(k,e) 6 E? : 0 < e _< £q and a e <_ k <_ k} 

•k 

where a > 0 is some constant such that 


(2.6) (k/e )X i (t ,e ) € ft for i = 1, • ♦ • ,m, OjCt_<T and (k,e) er'. 

Then there exist constants K > 0, 0<^5< 1 and kg <_ k, independent of t, 
k and e , such that all solutions of (1.8) satisfy the estimate 

(2.7) II y n H _< K5 ^ II II for all r-1 _< j < n < N and (k,e) £ T 
where r is the set defined on line (1.12). 


We make several remarks pertaining to the hypothesis in Theorem 2.1 before 
presenting its proof. 

1. If the eigenvalues of D^(t,e) lie in ft, that is, X^(t,e)€ ft for 
i = l,»»»,m, 0 <_ t _< T and 0 < e < £q, then there exist constants 

k, a and a trapezoidal region r so that condition (2.6) is 

satisfied. In this case let k be any positive constant and set 
* 

o=l. 

2. If condition (2.6) is satisfied for all 0 < e < Gq> then the set ft 
must contain the point at “. 

3. The condition that kg <_ k is natural in order to account for the 
rate of change of the eigenvectors of the matrix D^(t,e) which one 
misses by only examining the class of scalar test equations (1.1). 
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4. If the multistep method satisfies the conditions in Theorem 2.1 but is 
strictly n-stable and not strongly ft-stable, then estimate (2.7) 
does not hold. Approximate the equation y = — ( 1 / e ) y by the 

trapezoidal rule to justify this comment. 

The proof of Theorem 2.1 is a direct application of the following lemma 
which relates the behavior of a multistep approximation of a system of o.d.e.s 
with variable coefficients to the behavior of the multistep approximation for 
a class of scalar test problems (1.1). 


Lemma 2.2. Let T' denote the set defined on line (2.5). Consider the 
difference equation 

(2.8) y n = I Q i (k,e,x 0 (t n ),...,T r (t n 

j=l 

and (k,e) € T' 


))y n J for n = r,r+l, 


where y 11 is a m -component vector and each Qj, j = l,***,r, is an m x m 
matrix. Assume that the following two conditions are satisfied: 

1) There exists a constant Q, independent of t, k and e, such that 


(2.9) Qj(k,e,T 0 (t),**»,T r (t)) €C 1 (t,Q) for rk _< t < °° and (k,s)cr' 

and 


( 2 . 10 ) 




> T Q ( t ) , 1 


»T r (t)) = Qj(k,e,T 0 (t),«»*,x 0 (t)) 


kE. , 
J 


j = 1» 


where each matrix Ej, j = 1 , • • • , r , is uniformly bounded for all rk _< t _< » 
and all (k,e) £ r'. 
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2) For each point t satisfying rk < t < °° there exist constants 

*\j rv/ ^ 

K > 0 and 0 < 5 < 1, independent of t , k and e, so that the solutions of 
each "frozen coefficient" difference equation 


(2.11) v n = l 


l Q^k.e.x 
j = l J 


(t ), 


,e )> 




n = r , r+ 1 , ■ 


satisfy the estimate 


( 2 . 12 ) 


I v n II 


K6 (n_j) ll^ 


with 


~n r n n- 1 
v = (v ,v 


n-(r-l) ->T 
...,v ') 


for all r-1 < j < n < “ and (k,e) € T' . Then there exist constants 
K > 0, 0 ^ 6 < 1 and kg < 1c, independent of t, k and e, such that the 
solutions of (2.8) satisfy an estimate of the form 


(2.13) 


lly n || <_K6^ n ^ly J "l with y n = (y n ,y n 1 . 


n-(r-l) nT 

,y J 


for all r-1 < j < n < °° and (k,e) £ T where T is the set defined on line 

( 1 . 12 ). 


The proof of this result is lengthy and appears in Part II of the 
Appendix. 

Proof of Theorem 2.1 . The eigenvalues of the matrix (k/e^ 1 (x Q ( t) ,e) 
are distinct and are contained in fi for all 0 t _< T and (k,e) £ T' , and 
the raultistep method is strongly J2-stable, so 

bgk _ 1 

(2.14) ( a Q I + (— — )D^(ig(t) ,e) ) exists and is uniformly bounded 

for all rk < t < T and (k,e)£ T ' . 
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Consequently, for all (k,e) € T' the raultistep approximation of system (2.1) 
can be written as 


(2.15) 


r 



j=l 







where 


(2.16) 


R j (^ D ll (T 0 <t n ) ’ e) ’ (7) D u (T j (t n )>e)>k > £ ; 


b o\ \-i / b - k 


ro 1 + ( e ) D ll (T 0 (t n ),e) ) fcj 1 + ( e ) D 11 (l j ( t n ) ’ e) ) » 


j 


Simple calculations (see [27], Section 5 for details) establish that there 
exists a constant R, independent of k, t, and z, such that 


(2.17) R j ^(i)D 11 (T 0 (t),e), ] _ (x . (t) ,e)j £ C 1 (t ,R) 

and 

(2.18) 

= R j(^) D n (T o (t) « e) ‘ (7) D n (T o (t) ’ e) ) + 0(k) 

for all rk < t _< T and (k,e) e T' . 

•Jp 

Let t be any point in [ 0 , T ] and consider the "frozen coefficient" 


difference equation 
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(2.19) 


vI1= i ml R j(y D n (t > £) ’ y D n (t > £) j 




n = r, 


,N. 


/'k.'v * 

The eigenvalues of (^— j D ^ ^ ( t ,e) are distinct and contained in (2 for 
all (k,e) 6 T' and the multistep method is strongly , (2-stable, so the 
solutions of (2.19) satisfy an estimate of the form (2.12) for all t € [0,T], 

"k k 

Since 0 <_ Tg(t ) <^ T for rk <_ t <. T, the solutions of (2.11) satisfy 
an estimate of the form (2.12) for all t*€[rk,T], 

A direct application of Lemma 2.2 completes the proof of Theorem 2.1. 

Remark 2.1. If D^(t,e) €C^(t,e,B) for all 0 jC t < » and 0 _< e 
and satisfies the conditions in Assumption 1.2 for all 0 t < « and 
0 < £ _< Cq, and the matrix U(t,e) which diagonalizes D^(t,e) satisfies 


sup {llU(t,e)ll, HU 1 (t,e)n} < c 

0<t<- 
0 <£<£ 0 

where C is some constant, then estimate (2.7) holds for all j and n 
satisfying r-1 j n < °°. 

Remark 2.2 . A result in the spirit of this section was derived by Odeh 
and Liniger in Section III of [34], They consider one-leg approximations of 
nonlinear systems of o.d.e.s whch satisfy a monotonicity condition. They 
showed that if the one-leg method is A -stable (Odeh and Liniger [33]), then 

CO 

errors produced by initial data perturbations for the one-leg method decay 
exponentially with time. Analogous results were obtained for linear multistep 
methods by using the nonlinear transformation which relates the solution of a 
one-leg method to its corresponding linear multistep method. 
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The results of Liniger and Odeh cannot be used to obtain the result in 
this section because we consider systems (2.1) which do not necessarily 
satisfy the stated monotonicity condition and we consider a class of multistep 
approximations of systems (2.1) which is larger than the union of the one-leg 
and linear multistep methods. Even if we restrict our discussion to systems 
(2.1) which satisfy the stated monotonicity' condition for each fixed e 
satisfying 0 < e _< Eq and to the classes of one-leg and linear multistep 
methods, the results of Liniger and Odeh are still inapplicable in this 
section. The transformation which connects the solution of a linear multistep 
approximation of (2.1) to the corresponding one-leg approximation of (2.1) is 
unbounded as e 0 (van Veldhuizen [41]). 


3. PROPERTY 1.1 

In this section we apply the concept of stability domains together with 
the usual Dahlquist-Henrici stability theory to determine sufficient 
conditions under which a multistep approximation of (1.5) with the form (1.9) 
satisfies Property 1.1. 

Remark 3.1 . Throughout Section 3, Y will denote the trapezoidal region 
defined on line (1.12), we will assume that all constants are independent of 
the parameters k, t and e unless explicitly stated otherwise, and whenever 
we refer to systems (1.4) or (1.5) we will always assume that these systems 
satisfy the conditions in Assumptions 1.1 - 1.3. 

We begin by introducing various stability concepts for multistep 
approximations of systems of o.d.e.s which depend on a parameter. These 
definitions were previously stated and motivated in Section 3 of [27], 
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Cons ider the multistep method (1.9) and set v 11 = (v n ,**»,v n r+ ^)^ 
where v n is the solution of the corresponding homogeneous multistep method 
(1.8). The first definition corresponds to the usual meaning of stability in 
the sense of Dahlquist-Henrici. 

Definition 3.1 . 

a) The multistep method (1.9) is stable if for each fixed value of e 
satisfying 0 < e _< £q there exist constants a(e), K(e) and t(e) so that 
all solutions of the corresponding homogeneous multistep method (1.8) satisfy 
an estimate of the form 

a(e)t 

(3.1) llv n |l < K(e)e n ^ II II for all r-1 _< j _< n jC N 

when 0 _< k < t (e ) . 

b) The multistep method (1.9) is uniformly stable if it is stable with 
constants K(e), a(e) and r(e) which are independent of e for 

0 < e £q . 

Remark 3.2 . If the matrix A(t,e) is continuous for all 0 _< t _< T and 

0 < e £q, then method (1.9) is stable provided that the roots of the 

r 

characteristic polynomial p(X) = £ a.X^ satisfy the root condition. 

j=0 J 

Furthermore, if A(t,e) is continuous for all 0 ^ t _< T and 0 e _< e^, 
then this same algebraic condition implies that (1.9) is uniformly stable. 

The next definition introduces stability concepts which are appropriate 
for multistep approximations of stiff systems of o.d.e.s which depend on a 
parameter when the step size resolves only the slow time scales in the system 
of o.d.e.s and not the fast time scales. 
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Deflnition 3.2. 

a) The multistep method (1.9) is uniformly (k,e)-bounded in r if 

(a^I + bgkA(xQ(t) ,e)) 1 exists and is uniformly bounded for all 

(3.2) rkj(t^T and (k,e) £ F, 

and there exists a constant K so that all solutions of the corresponding 
homogeneous raultistep method (1.8) satisfy an estimate of the form 

(3.3) II v n || _< Kllv 11 1 II for all r _< n _< N and (k,e) € T. 

b) The multistep method (1.9) is uniformly (k,e)-stable in T if 

condition (3.2) is satisfied and there exist constants K and a so that all 
solutions of the corresponding homogeneous multistep method (1.8) satisfy an 
estimate of the form 

(3.4) nv n || £ Ke n J II v~* || for all r-1 £ j _< n _< N and (k,e) £ r. 

c) The multistep method (1.9) is strongly (k,e)~stable in r if 

condition (3.2) is satisfied and there exist constants K and 0 < 6 < 1 so 
that all solutions of the corresponding homogeneous multistep method (1.8) 
satisfy an estimate of the form 


(3.5) ||v n || < K5^ n ^llv^n for all r-1 < j < n < N and (k,e) € T 
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The following lemma states an obvious relationship among the concepts 
defined in Definition 3.2. 

Lemma 3.1. 

a) If (1.9) is strongly (k,e) ~stable in r, then (1.9) is uniformly 
(k,e) -stable in T . 

b) If (1.9) is uniformly (k,e) -stable in r, then (1.9) is uniformly 
(k,s ) -bounded in r. 

Remark 3.3. 

1) The homogeneous multistep method (1.8) satisfies Property 1.1 if and 
only if there exists a trapezoidal region r with the form (1.12) such that 
(1.8) is uniformly (k,e)-stable in r. 

2) In Theorem 2.1 we used the geometry and properties of stability 
domains to determine sufficient conditions under which there exists a 
trapezoidal region with the form (1.12) such that a multistep approximation of 
system (2.1) is strongly (k,e)-stable in r. 

In order to precisely state the main results which appear in Theorems 3.1, 
3.2, and Corollary 3.1, we find it convenient to separate implicit homogeneous 
multistep methods (1.8) into three different classes. 

Class X ; b^ = 0, i = l,»**,r. 

Class Y : -r^(t) = Tp(t) for i = l,«««,r, and 0 _£ t _< T. 

Class Z : The integers i = l,«««,r, can be partitioned into two sets 

1^ and I 2 such that 

i £ 1^ if b^ = 0 or x i (t) = Xq( 0 f° r a H c e [ 0 , T ] 


(3.6) 
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and 

(3.7) i € if * 0 and T^(t) t tg(t) for all t£[0,T]. 

It is easy to check that Class X consists of the backward difference 
methods. Class Y consists of the one-leg methods and Class Z consists of 
multistep methods which are neither backward difference methods nor one-leg 
methods. 

Let T(t,e) denote the matrix which appears on line (1.6). For 
i = 1 , • • • , r, let E^(k,t,e), E* 2 (k,t,e), E 21 (k,t,e), E 22 (k,t,e) € C p ^t.e.B) 
denote matrices with dimensions m x m, m x (n-ra), (n-m) x m and 
(n-m) x (n-m), respectively. 

Theorem 3.1. (Property 1.1 for homogeneous multistep methods (1.8).) 
Approximate system (1.4) by a multistep method (1.8). Assume that the 
multistep method is implicit, stable, the step size k resolves the slow time 
scales in (1.4), that is k _< kg where kg is any constant satisfying 

(3.8) k 0 1 b 0 1 ^ * for 0 < t < T and 0 S. e S. e 0’ 

and the step size does not resolve the fast time scales in (1.4) but the fast 
time scales are attenuated, that is 


method (1.8) is a strongly (k,e) ~stable approximation of the equation 
dw ^ 

(3.9) e dt" = D^(t,e)w _in_ r where a is some constant and kg is 
defined by (3.8). 
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If (1.8) lies in Class X or Y, then (1.8) is uniformly (k,e)-stable in 

r. 

If (1.8) lies in Class Z and the set I 2 is empty, then (1.8) is 

uniformly (k,e)-stable in r. 

If (1.8) lies in Class Z, the set I 2 is not empty and for each ie^ 
the matrix T(t,e) satisfies the coupling condition 

(3.10) T - 1 (T 0 (t),e)X(T i (t), e ) 
and 

(3.11) (k, t , e) = 0(e) 

for all rk £ t <_ T and 0 <_ e €q, then (1.8) is uniformly (k,e)-stable in 

r. 

If (1.8) lies in Class Z and there exists an integer i € I 2 so that 

condition (3.10) is satisfied but condition (3.11) is violated, then (1.8) 

need not be uniformly (k,e)-stable in r for any choice of the constants 
* 

a and kg . 

The central idea behind the proof of Theorem 3.1 is to separate the scales 
of the multistep approximation of (1.4). We show that after an appropriate 
bounded change of variables, a given multistep approximation of (1.4) can be 
transformed into the same multistep approximation of the system = D(t,e)y 
plus a perturbation. (D(t,e) is the block diagonal matrix which appears in 
Assumption 1.1.) If the hypotheses in the first paragraph of Theorem 3.1 are 
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satisfied, then we can uniformly bound the multistep approximation of 

= D(t,e)y for all (k,e) € T. We use stability to bound the solution of 

dy II 

the multistep approximation of the system — = D 22 (t ,e)y-j- ^ and we apply 

assumption (3.9) to bound the solution of the multistep approximation of the 
dy I 

system e = D^(t,e)y^.. Furthermore, in the first four cases stated in 
Theorem 3.1 we can uniformly bound the perturbation and apply the following 
variant of a theorem due to Strang [38] to estimate the solution of the 
difference method. 


Lemma 3.2. Consider the difference approximations 


(3.12) 


y 11 = l Q j (x 0 (t n ),...,T r (t n ),k, i 

j=i 


)y n J , n = r , • • • ,N, 


and 


(3.13) 


r „ 


w n = l Q.(x 0 (t n ), ...,x r (t n ),k,e)w n J , n = r,.**,N, 

j=l 


where Qj, Qj, j = 1 , • • • , r , are n x n matrices . If the difference 

approximation (3.12) is uniformly (k,e)~stable in r and 


(3.14) 


Qj = Qj + 0(k) for all (k,e) € r and t r <. t n <. T, 


then the difference approximation (3.13) is also uniformly (k,e) -stable in 

r. 


Proof of Theorem 3.1 . The matrix ( a Q^ + bQkA(xQ( t ) , e) ) ^ exists and is 
uniformly bounded for all rk < t _< T and (k,e) € T, so (1.8) can be written 


in the form 
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(3.15) 

where 


r 

v 11 = l R.(A(T 0 (t n ),e),A(T.(t n ),e),k,e)v n 2 , n = r,...,N, 

j=l J J 



j = 1, 




A simple application of identity (1.6) shows that 


Rj ( A(xg(t) ,e ) >A(t j (t),e),k,e) 


= T(x 0 (t),e)A 1 T 1 (x 0 (t),e)T(x j (t),e)A 2 T ^x^O.e) 


where 



" ) 

(^qI + ^QkD22( T o(t) ,e ) ) ! 


and 



Make the change of variables v 11 = T(x n (t ),e)w n in (3.15) to obtain the 

U n 


difference equation 
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(3.16) 
where 

(3.17) 


w n = Y R. (k,t ,e)w n ^ 
J n 

J = 1 


Rj (k, t ,e ) = T _1 (To(t),e)Rj (A(x 0 (t) ,e) ,A (tj (t) ,e) ,k, e)T(-r 0 (t-jk) , e) , 


j = !»•••, r. 


Let D(t,e) be the matrix defined on line (1.6) and set 
Rj{D(x 0 (t) , e),D(tj ( t) ,e),k,e) 

(3.18) 

= (a Q I + b Q kD(T 0 (t) ,e)) 1 (a^ I + kb^ D(ij (t ) ,e) ) , j = l,***,n. 
If method (1.8) lies in Class X, then 
Rj (A.(TQ(t) ,e ) ,A(Xj (t) ,e) ,k, e) 

= (a Q I + b Q kA(T 0 (t) ,e)) ^a^I), J = l>***>r, 

and 


(3.19) 


Rj(k,t,e) = Rj (D(x 0 (t) ,e) ,D(tj (t) ,e) ,k,e) + 0(k) 

for j = !,•••, r, rk < t < T and (k,e)er. 


If method (1.8) lies in Class Y, then 
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Rj(A(T 0 (t),e),A(Tj(t),e),k,e) 


,A ii 0 


= T(T Q (t),e) 


|T l (T 0 (t),e), j - !»•••, r, 


0 A 


22 > 


where 


A^Ct.e) 


= (v + HH D ii (T o (t) ’ e) ) -1 fj 1 + ^*T") D ii (T j (t) > e) )’ 


j = !,•••, r. 


and 


A 22 (t,e) = (a Q I + b 0 kD 22 (T 0 (t),e)) 1 ( a j I + b j kD 22 ( T j (t) ,e) ) , 


j 


= 1 , 


r. 


a 

Consequently, Rj(k,t,e) satisfies (3.19). 

If method (1.8) lies in Class Z and the set I 2 is empty, then 
Rj(k,t,e) satisfies (3.19). 

A 

If method (1.8) lies in Class Z and j £1^, then Rj(k,t,e) satisfies 
(3.19) for all j £ 1 ^ and (k,e) £ T. If j £ I 2 and condition (3.10) is 

satisfied, then 


where 


R (A(T 0 (t),e),A(T j (t),e),k,e) = y| + Y 2 + Y 3 
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and 


with 


A 


T(T 0 (t),e) 


|T _1 (t j ( t ) , e) , 


' 22 ' 


Y 2 = kT(T Q (t),e)A 1 | 


r E J u (k,t,e) Ej 2 (k, t,e)^ 


E^ 2 (k,t,e) y 


U 2 T V^O.e), 


0 O' 


Y3 = kT(T Q (t),e) j 


| T _1 ( T j (t) ,e) 


*21 °- 


-1 • / kb. * 

3 2 1 " (sq 1 + kb 0 D 22 (T 0 (t),e)) E^k.t.e) (a^I + (-^-) D 1 i (tj ( t) , ej 


and y 2 are uniformly bounded for all (k,e) € T. y 2 is uniformly 
bounded if condition (3.11) is satisfied. Consequently, Rj(k,t,e) satisfies 
(3.19) when (1.8) lies in Class Z and the coupling condition (3.10), (3.11) is 
satisfied for all i£I 2 . 

By Remark 3.2, method (1.8) is a uniformly stable approximation of the 
equation = D 22 (t,e)y for all 0 <_ t _< T and (k,e) eT. This fact and 
assumption (3.9) together imply that the difference equation 


(3.20) x n = l R.(D(T n (t),e),D(T.(t),e),k,e)x n ^ , n = r,*«»,N, 
j=l J J 


is uniformly (k,e)-stable in r 
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Identity (3.19) holds for the first four cases in this theorem, so we can 
conclude that the difference equation (3.16) is uniformly (k,e)-stable in F 
by Lemma 3.2. 

The matrices T(xg(t),e) and T *(xg(t),e) are uniformly bounded, so we 
can conclude that (3.15) is uniformly (k,e)-stable in r and the proof of 

Theorem 3.1 in the first four cases is complete. We establish the last 
statement in Theorem 3.1 by constructing an explicit example. This example is 
presented in Section 5. 

Theorem 3.2. (Property 1.1 for inhomogeneous multistep methods (1.9).) 
Approximate system (1.5) by a multistep method (1.9) and set 

u n = (y n ,***,y n r+ ^)^. Assume that method (1.9) is implicit, stable, the step 
size k satisfies 0 < k kg where kg is defined on line (3.8), and 
(1.9) satisfies assumption (3.9). 


If the corresponding homogeneous difference method (1.8) lies in Class X 
or Y and F(xg(t ) , • • • ,t ( t ) ,k, e) = F(xg(t),e), then there exist constants K 
and a such that 


(3.21) 


w i- . . 

Ilu n || < Ke n 3 fllu 3 ll + t .f) for all 
— n-j 1 


r-i < j < n < N 
and (k,e) er 


where 

f = sup ll(f T (t,e),f TT (t,e)) T H. 
0<t<T 


0<e<e 0 
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That is, method (1.9) satisfies Property 1.1. 

If the corresponding homogeneous difference method (1.8) lies in Class Z, 
the set I 2 is empty and F(xg(t) , • • • ,x r (t) ,k,e) = F(xg(t),£), then all 
solutions of (1.9) satisfy an estimate of the form (3.21). 

If the corresponding homogeneous difference method (1.8) lies in Class Z, 
the set I 2 is not empty, conditions (3.10) and (3.11) are satisfied for 
each i € I 2 and F(xg(t) , • • • ,x r (t) ,k, e) is a linear combination of 
F(x.(t),e) with i £ I 2 , then all solutions of (1.9) satisfy an estimate of 
the form (3.21). 

Assume that the corresponding homogeneous difference method (1.8) lies in 
Class X, Y or Z, conditions (3.10) and (3.11) are satisfied for each i € I 2 
if (1.8) belongs to Class Z, and that F(xg(t) , • • • ,x (t) ,k,£) is a linear 
combination of F(x^(t),e), i = 0,**»,r. If conditions (3.10) and (3.11) are 
also satisfied for all i£I^, then all solutions of (1.9) satisfy an estimate 
of the form (3.21). If condition (3.10) is also satisfied for all i £ Ij but 
condition (3.11) is violated for some value of i £1^, then the solutions of 
(1.9) need not satisfy an estimate of the form (3.21). Consequently, the 
inhomogeneous multistep method (1.9) need not satisfy Property 1.1 in this 
last case. 

The proof of Theorem 3.2 depends on the discrete version of Duhamel's 
Principle and the following well-known lemma: 

Lemma 3.3. For any values of k and e such that 

(a^I + kbgA(xg(t) ,e) ) * exists for 0 _< t £ T, 


the multistep method (1.9) can be written as 
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(3.22) 


V n = l R,(A(T 0 (t n ),e),A(T ,(t n ),e),k,e)y n J 

j=l J J 


+ kF(T 0 (t n ),.*. > T r (O,k,e), n = 


where 


(3.23) 


R. = (a Q I + kb 0 A(T 0 (t n ),e)) ^a.I + kb .A(t .(t n ), e )) 


j = 1 , • • • »r , 


and 

(3.24) F = (a Q I + kb Q A(T 0 (t n ),e)) 1 F(t Q (t n ) , . • • ,x r (t n ) ,k,e) . 


Furthermore, (3.22) can be written as the equivalent one-step method 


(3.25) 


u n = a(k,t ri> e)u 11 1 + kG(k,t n ,e), 


n = r. 


,N, 


where 


(3.26) 


R (a(t (t ) ,e ) , A (t (t ),e),k,el***R (A(x (t ),e),A(x (t )e),k,e 
1 0 n In r On rn 


a = 


I 

0 


0 • • • 

I • • • 


0 

0 


0 • • • •! 


I denotes the n x n identity matrix and 
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(3.27) u n = (y r 


,y n r+1 ) T ,G(k,t n ,e) = F('r 0 (t n ),***,T r (t n ),k, e ),0,«.. > 0) T . 


Proof of Theorem 3.2. We first note that in all cases listed in this 
theorem, method (1.9) is uniformly (k,e)-stable in r by Theorem 3.1. 

Method (1.9) can be written as the equivalent one-step method (3.25) - 

(3.27) . For any integers n and j satisfying r-1 _< j _< n _< N, the 

discrete version of Duhamel's Principle implies that the solution of this one- 
step method is given by 

n 

(3.28) u n = S(n,j)u J + k £ S(n, j)G(k,t . ,e) 

i-j+1 J 

where 

n 

/ II a(k,t.,e) if n _> j + 1 

J i-j+1 1 

S(n, j ) = \ 

( I if n = j 


Estimate (3.28) with the triangle inequality to obtain 

n 

(3.29) || u n II _< || S ( n , j ) I! llu J II + k £ BS(n,j)|| llG(k, t . ,e) II . 

i-j+1 J 

In all cases listed in this theorem (1.9) is uniformly (k,e)-stable in 
T, so there exist constants and a so that 

at 

(3.30) II S ( n , j)ll e n ' 1 for all r-1 _<j an< i (k,e) € r. 
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Consequently , the solution of (1.9) will satisfy an estimate of the form 
(3.21) if there exists a constant C 2 such that 

(3.31) nG(k,t,E)n < C 2 f for all 0 _< t _< T and (k, £ )£T. 

By identities (3.24) and (3.27), inequality (3.31) holds if 

(3.32) Hp|| = || ( a Q I + b 0 kA(x 0 (t),e)) 1 F(x Q (t ) , • • • ,x r (t) ,k,e) II < C 2 ? . 

Set 

Mj = (a Q I + (b 0 k/e)D 11 (x 0 (t),e)) 1 and M 2 = (a Q I + b Q kD 22 (x 0 (t), 0 ) ) *. 

If (1.9) satisfies the conditions in any of the first three cases listed in 
this theorem, then 


- 1 , 


F = (a Q I + b Q kA(x o (t),0)) F(xq( 0 ,e) 


T(x 0 (t),e)| 


0 


0 M 2, 


|T" 1 (x o (t),0)F(x o (t), e ) 


T(xg(t) ,e) 


/ M 1 

°v 


\ 0 

hJ\ 
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M 1 

The matrices ( — ), and T(xQ(t),e) are all uniformly bounded, so we can 

conclude that f satisfies an estimate of the form (3.32). 

If (1.9) satisfies the conditions in the remaining cases listed in the 
theorem, then 


F = (a Q I + b Q kA(T 0 (t) ,e)) 1 / £ cuFtx^O.e) 

\ i=0 i 


v-1. 


= l a jL (a () I + b 0 kA(x 0 (t),e)) F(T i (t),e) 
i=0 


where cu, i = (),•••, r, are constants. For each value of i, 


-1 


(a Q I + b Q kA(T 0 (t) ,e)) F(r ;L (t),e) 


= T(T 0 (t),e)| 


'Mj 0 


0 M, 


|T 1 (T 0 (t),e)T(x i (t),e)T 1 (x i ( t ) , e)F(x i (t) , e) 


M, 


(“ X 1 + kE^) kM^J 


12 


fl(Xi(t),e)' 


= T(x 0 (t),e) 


■«2(-r) 


M 2 (I + kE * 2 )/ 


All matrices which appear on the previous line are uniformly bounded if 

condition (3.11) is satisfied. If condition (3.11) is violated, then 
E"*" 

2 1 1c 

kM„( ) = 0(— ) in general and the proof of Theorem 3.2 is complete. 

“ £ £ 
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Remark 3.4. 

1) In [29] we construct an explicit example of a multistep method which 
satisfies all conditions in the last statement of Theorem 3.2 and does not 
satisfy Property 1.1 for a particular system of o.d.e.s satisfying Assumptions 
1.1 - 1.3. Consequently, the results of Theorem 3.2 are sharp. 

2) Strictly speaking, an r-step multistep approximation of system (1.5) 
satisfies Property 1.1 if: 

i) The solutions of the multistep method satisfy an estimate of the form 
(3.21). 

1 r -l 

ii) The (r-l) additional values of the r-step method, y ,•••, y , 

are computed by a one-step method whose solutions satisfy an estimate of the 

form (3.21) for 0 <_ j n £ r-l (with u^ replaced by yJ). 

Consequently, the additional initial values of a multistep method must be 
appropriately computed in order for a multistep method to satisfy Property 
1.1. This point is discussed on pp. 39-40 in [29]. 

In order to compare the properties of one-leg and the classical linear 

multistep methods, we restate Theorems 3.1 and 3.2 in terms of these methods. 

Corollary 3.1. Approximate system (1.5) by a multistep method (1.9). 
Assume that method (1.9) is implicit, stable, the step size k satisfies 

k <_ k^ where kg is defined on line (3.8), and (1.9) satisfies assumption 
(3.9). 

If the multistep method is a one-leg method or a backward difference 

method, then it is uniformly (k,e)-stable in T and all solutions of the 

inhomogeneous one-leg or backward difference method satisfy an estimate of the 

form (3.21). 
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If the multistep method is a linear multistep method which is not a 
backward difference method and the matrix T(t,e) which appears on line (1.6) 
satisfies the coupling condition 


(3.33) T 1 (t,e)T(t-k,e) = I + k 


^ E 1 i (k> t > g) E 1 12 (k,t,e)^ 


k E 21 (k,t,e) E*,(k,t,e) y 


22 ' 


and 


(3.34) eJj = 0(e) 

for all 0 < t < T and 0 < e < Cr, where e!., i,i = 1,2, are characterized 

— — — — 0 l j ’ 

in Theorem 3.1, then the multistep method is uniformly (k,e)-stable in r. 

Furthermore, all solutions of the inhomogeneous linear multistep method 

satisfy an estimate of the form (3.21). 

If the multistep method is a linear multistep method which is not a 

backward difference method and condition (3.33) is satisfied but condition 

(3.34) is violated, then the linear multistep method need not be uniformly 

* 

(k,e)-stable in T for any choice of the constants a and kQ. 

We now discuss the coupling condition (3.33), (3.34) before making some 

remarks about the practical implications of Corollary 3.1. 

Assume that the matrix T(t,e) which appears on line (1.6) satisfies 
2 

T(t,e)eC (t,e,B). Make the change of variables y = T(t,e)w in system (1.5) 
to obtain the equation 

= ^D(t,e) - T _1 (t,e) ' — < - ^ ) w + T _1 (t,e)F(t,e) 


(3.35) 
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where the matrix D(t,e) is defined on line (1.6). Identity (3.35) shows 

that the matrix - |t ^(t,e) controls the coupling between the fast 

and slow time scales in system (1.5). 

Now consider the expression on line (3.33). By Taylor's Theorem 


T(t-k,e) = T(t,e) - k dT( ^’ e) + 0(k 2 ) 

(3.36) 


for 0 t _< T and 0 _< e e 0* 


Substitute (3.36) into (3.33) and let k + 0 to obtain 


(3.37) 


T~ ^ ft- r 1 ^ ^ e ) — 

-T (t,e) ^ 


'E^O.t.e) E} 2 (0,t,e)' 


E 21 ( 0 ,t,e) £22(0,1,6), 


If condition (3.34) is also satisfied, then 


(3.38) -T~ X (t ,e ) dT( ^’ e) - - 


^(O.t.e) E} 2 (0,t,e) 


E 22 (0, t , e) 


+ 0(e). 


Identity (3.38) shows that conditions (3.33) and (3.34) impose a 
restriction on the strength of the coupling from the fast to the slow time 
scales of system (1.5). In particular, conditions (3.33) and (3.34) are 
satisfied if the coupling from the fast time scales of (1.5) to the slow time 
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scales of (1.5) is sufficiently weak. For this reason conditions (3.33) and 
(3.34) are called the coupling condition in this paper and in [27], A 
detailed discussion of the coupling condition can be found in van Veldhuizen 
[41]. 

Condition (3.38) clearly places a restriction on the structure of system 
(1.4) and the matrix T(t,e) in particular. If (3.38) is satisfied, then 
standard results from the theory of o.d.e.s imply that 


(3.39) 


M u (t,e) M 12 (t,e) 

T(t,e) = T(0,e) I I + 0(e) 

\ n \ , 

22 ' 


0 M (t,e) 


for 0 t T and 0 S. e e Q 


where Mjj is an m x m matrix, M [2 is an m x (n-m) matrix, and M 22 is 
an (n-m) x (n-m) matrix. We obtain the following result: 


Lemma 3.4. _If_ 

1) the matrix A(t,e) in (1.4) is constant 
or 

2) system (1.4) appears in singular perturbation form and has no turning 
points, that is, 


'A^tje) A 12 (t,e) 


A(t,e) = 


A 21 (t,e) 


A 22 (t,e) 
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where A.^t.e) eC P (t,£,B) for i, j = 1,2, and A^Ct.e) is invertible with 
A^(t,e) 6C^(t,e,B), then condition (3.39) is always satisfied . 

The proof of Part 1) of Lemma 3.1 is obvious and a proof of Part 2) can be 
found in Section 2 of [28], 

It is well-known that a one-leg method and its corresponding linear 
multistep method have identical properties when they approximate systems of 

0. d.e.s with constant coefficients. Furthermore, there is an invertible 
transformation which transforms the solution of a linear multistep method into 
its corresponding one-leg twin, Dalhquist [8]. Despite these facts, Dahlquist 
and Lindberg [11], Dahlquist [10], Nevanlinna and Liniger [32], and van 
Veldhuizen [41] have demonstrated that in various ways the behavior of a 
linear multistep method and its corresponding one-leg method can be very 
different when they approximate variable coefficient (and nonlinear) stiff 
problems. We make the following additional observations about these two 
classes of methods on the basis of Corollary 3.1: 

1 . The uniform (k ,e )-stability (stability with large step sizes) of linear 
multistep approximations of stiff linear systems depends on the form in 
which system (1.5) appears . This comment is not true for one-leg 
methods. Consider the cases when the matrix A(t,e) in (1.5) is constant 
or system (1.5) appears in singular perturbation form without turning 
points, that is, system (1.5) has the form 

dy I 

e = A u (t,e)y I + A 12 (t,e)y II + f ] .(t,e) 

= A 21 (t,e)y I + A 22 (t,e)y II + f Ii; (t,e) 


dy II 

dt 
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where the matrices A_(t,e), i = 1,2, j = 1,2, are characterized in Lemma 
3.4 and f^(t,e), f^^Ct.e) € C p (t,e,B) . In these two cases the uniform 
(k,e )-stability of both a one-leg approximation of system (1.5) and its 
corresponding linear multistep approximation are determined by the usual 
Dahlquist-Henrici stabilty theory to control the slow time scales of (1.5) 
and the use of stability domains to control the fast time scales. If 
system (1.5) satisfies the conditions in Assumptions 1.1 - 1.3 but is 
otherwise unrestricted, then the two previously mentioned properties 
determine the uniform (k, e )-stability of a one-leg approximation of 

(1.5) . The uniform (k,e)-stability of a linear raultistep approximation 

of system (1.5) depends on an additional condition, namely the coupling 
condition (3.33), (3.34). This condition only arises when the 

coefficients in system (1.5) vary with time. 

We note that if we make the change of variables y(t) = T(t,e)w in 

(1.5) , then the resulting equation for w appears in singular perturba- 
tion form. According to Corollary 3.1, the uniform (k, e)-stability 

properties of a linear multistep approximation of the scaled system 

which w satisfies are better than the uniform (k,e)-stability 

properties of the linear multistep approximation of the original system 

(1.5) . 

2. If we consider problems which exhibit a strong coupling from the fast time 
scales to the slow time scales, then a one-leg method using a given step 
size k can compute an accurate solution of the equation while the 

solution of its corresponding linear multistep method using the same step 
size will produce a solution which is not accurate at all. In other 
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words, the linear multistep method requires a smaller step size to obtain 
the same accuracy as the corresponding one-leg method. This point is 
illustrated with an example in Section 5. 

Remark (2) is consistent with observations made in [14] and [15] which 
demonstrate that the integration package TRAPEX (based on the trapezoidal 
rule with extrapolation) has difficulties solving nonlinear systems of 
o.d.e.s when there is a strong coupling between the slow and fast time 
scales. According to Corollary 3.1, these numerical difficulties could be 
avoided if the trapezoidal rule is replaced by the implicit midpoint 
rule. This change was made by Lindberg [25], [26] in the integration 

package IMPEX. 

3. Corollary 3.1 demonstrates that one-leg methods have better stability 
properties with large step sizes, in general, than their linear multistep 
counterparts. Consequently, one-leg methods are to be preferred over 
their linear raultistep counterparts when choosing a discretization for a 
system of stiff o.d.e.s. 


4. THE RELATIONSHIP BETWEEN PROPERTY 1.1 AND D-STABILITY 

In this section we relate the theory developed in Section 3 to the concept 
of D-stability introduced by van Veldhuizen in [41], In [41] van Veldhuizen 
stated that D-stability complements the usual notions of stability and 
stability domains, and that D-stability has no relation to the propagation of 
errors for solutions of difference methods. An important consequence of the 
results in this section is that we relate the concept of D-stability to the 
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usual notions of stability and stability domains, and to the propagation of 
errors for multistep methods which use large step sizes. 

Following van Veldhuizen [41], consider the homogeneous model problem 

(4.1) = A(t)y , t > t Q , 

where A(t) is a. matrix. Approximate system (4.1) by a discretization method 
which can be written in the form 

(4.2) y ±+1 = G(t i ,k)y i 

where k is the step size, = t Q +ik, and G(t is k) is a matrix. If the 
discretization is a one-step method, then y^ is an approximation to 
y(t^). If the discretization is a multistep method, then the recursion (4.2) 
describes the canonical one-step recursion corresponding to the multistep 
method . 

Let T, kg and X denote given constants. 

Definition 4.1 . (van Veldhuizen) The discretization of (4.1) resulting 
in the recursion (4.2) is called D(£)-stable if for all stiff systems in the 
Class £, for all t g [tg,tg+T] and all k£(0,kg] 

II G( t ,k) II <_ M < °° 

where M is a constant depending only on kg and the Class Z. 


Van Veldhuizen makes the following choice for the Class Z: 
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Definiton 4.2 . (Class £ ) The Class £ of stiff systems consists of all 
linear systems (4.1) parametrized by a parameter e 6 (0,Eq] which satisfy the 
following conditions: 

(51) y(t) € (j) 2 

(52) A(t) = T(t)AT _1 (t) for all te[t 0 ,t 0 +T]. Here A(t) is a 
diagonal matrix given by 


A(t) 




L o 


o 

n(t)_ 


where Re(A(t)) _< A < 0 for all t € [tQ,tg+T]. 

(S3) A, p, T, T ^ depend smoothly on t and possibly e €(0,Eg], and 

the derivatives from order zero up to a certain order (suf f iciently high) are 

bounded uniformly in t and e C(0,Eq]. 

Discretizations of systems (4.1) fall into two categories. Either all 

time scales in (4.1) are resolved by the step size or some time scales in 

(4.1) are not resolved by the step size. (For systems which lie in Class £, 

all scales are resolved if — + 0 as k + 0. The fast time scales are not 

e 

k * * 

resolved if — > a with a a constant as k + 0.) If all scales in (4.1) 
£ — 

are resolved by the step size, then the usual Dahlquist-Henrici stability 
theory controls the error propagation properties of the multistep 
discretization. A meaningful stability concept for a multistep discretization 
of stiff equations which complements the Dahlquist-Henrici theory should apply 
to the case when some time scales in (4.1) are not resolved by the step 
size. Consequently, we alter van Veldhuizen's Definition 4.1 for systems 
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(4.1) which depend on a parameter 0 < e < to reflect the complementary 

case when all scales in (4.1) are not resolved by the step size. 

Let r denote the trapezoidal region defined on line (1.12). 

Definition 4.1a . Consider systems (1.4). The discretization of (1.4) 
resulting in the recursion 

(4.2a) y i+1 = G(t i ,k,e)y 1 

is called D(£)-stable if for all stiff systems in the Class £, for all 
t € [ t Q , t Q +T ] and (k,e)€T 


llG(t,k,e)H £ M < ~ 

■k 

where M is a constant depending only on and the Class £. 

We also replace the Class JC defined in Definition 4.2 by the Class C 
which consists of all systems (1.4) (or a subset of systems (1.4)) which 
satisfy the conditions in Assumptions 1.1 and 1.2. 

We are now prepared to state two results which connect our version of van 
Veldhuizen's Definition 4.1 to the results in Section 3 and to the propagation 
of errors for multistep methods which use large step sizes. 

Theorem 4.1. Approximate a system in Class C by a multistep method 
(1.8). The multistep method is D(C) -stable if and only if for all stiff 
systems in C, the multistep method is uniformly (k, e)~bounded in T with a 


bound which only depends on the Class C 
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Theorem 4.2. If the hypotheses in Theorem 3.1 hold uniformly for all 
systems in Class C, then the multistep method is D(C) -stable. Furthermore, 
for all systems In C the multistep method is uniformly (k > e) ~stable In r 
with constants K and a (which appear on line (3.4)) depending only on the 
Class C. 


We have applied Part b) of Lemma 3.1 to justify the first statement in 
Theorem 4.2. 


5. NUMERICAL EXAMPLES AND THE PRACTICAL IMPLICATIONS OF PROPERTY 1.1 

In this section we introduce a model variable coefficient system of 
o.d.e.s which satisfies the conditions in Assumptions 1.1 - 1.3. We 

approximate this model equation by a particular one-leg method and its 
corresponding linear multistep method. Our goal is to illustrate the 
practical consequences of Corollary 3.1 which are stated at the end of Section 
3, and to make a general comparison between the classes of one-leg and linear 
raultistep methods. 

Consider the model system 


(5.1) 


where 


dv 

dt 


o 

rr 

<», 0 < e < • 1 , 

v(0) = ( 

cos(t) 

sin(t) \ i 

H °\ 




sin(t) 

cos(t)/ 1 

( 0 -l/ 


(5.2) 


U(t) = 
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and U (t) denotes the transpose of U(t). 

System (5.1) satisfies the conditions in Assumptions 1.1 and 1.2 with the 
matrices T(t,e) and D(t,e) which appear in Assumption 1.1 given by 


(5.3) 


T(t,e) = U*(t) 


and D(t,e) = D. 


The matrix U(t) defined on line (5.2) is unitary, so 


T 1 (t ,e )T(t-k,e ) = U(t)U*(t-k) = 


cos(k) sin(k)' 


,-sin(k) cos(k) 


= I + k 


0 1 


(-1 0 , 


+ 0 ( 0 . 


Consequently, system (5.1) does not satisfy the coupling condition (3.33), 
(3.34). 

System (5.1) can be transformed into a system of o.d.e.s with constant 
coefficients. Make the change of variables v = U (t)w in (5.1) to obtain 


(5.4) 



-1 


w, w(0) = v(0), 


-1 


Equation (5.4) is easily solved to produce the general solution of (5.1) which 
is given by 


v(t) = U (t)M 




M *v(0) 
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where 


2e 

1 + J - e 


M = 


2 e 


» X ± = V 2 (“(I + 7 ) ± 7 J ) • and J = 


2e - 3e 


1 + J - e 


A simple calculation shows that 


X +0(1) and X , 

e + 


= -1 + 0(e) as e 0. 


We conclude that the solutions of (5.1) behave like the solutions of the 
constant coefficient system = Dv. Furthermore, any solution of (5.1) can 
be written as the sum of a slowly varying solution of (5.1), denoted by 
v s (t), and a rapidly varying solution of (5.1), denoted by v^(t), namely, 


v(t) = v S (t) + V R (t) 


where 


v“(t) = U (t)M 


M ^v(O) and v R (t) = U (t)M 


0 e 


x+t 


fe X - t 0 


0 0 


M ^v(O) , 


Consider the one-leg approximation of (5.1) given by 


(5.5) 


i 1 - kU *(! c „ t l + 7 'nWl t n+l + 7 t j ) x ” +1 

* (1 + *“*(7 t n+l + 7 t „)' ,u (7 t „+i + 7 '„))*”> *° * v <°>' 
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The corresponding one-step approximation of (5.1) is given by 

(5.6) (I - U*(t n+1 )DU(t n+1 ))y n+1 = (I + j U*(t n )DU(t n ))y n , y° = v(0). 

Let T denote the trapezoidal region 

(5.7) T = { (k,e) : 0 < e < .1 and a*e < k < k Q } 

* 

where a is any positive constant and kg is any constant satisfying 
kg < 4/3. Both the one-leg method (5.5) and its corresponding one-step method 
(5.6) are implicit, stable, and strongly (k,e)-stable approximations of 
s y = -y in T. Corollary 3.1 implies that (5.5) is uniformly (k,e)-stable 
in T . 

We now show that method (5.6) is not a uniformly (k,e)-stable 

* 

approximation . of system (5.1) in T for any choice of the constant a . 
Since method (5.6) satisfies all hypotheses in the first paragraph of 
Corollary 3.1 but system (5.1) does not satisfy the coupling condition (3.33), 
(3.34), we conclude that the results in Corollary 3.1 (and Theorem 3.1) are 
sharp. 

The matrix (i - U (t)DU(t)) is invertible for all t 0 and 
(k,e) er, so method (5.6) can be written as 
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; r, + l> DU( Vl 


» ‘(I + T 0*(t )DO(t ))y” 


■ "X+iXi - r- + ! D)u<t n )y" 

(5.8) 



where and ^22 are un if° rm ly bounded in r and 

M 21 = ^ - 4|) cos(t n + l )cos(t n ) ^ sin(t n )c0s(t n+l ) 

- sin(t . ) cos ( t )1 + E 
n+1 n / 


where E is uniformly bounded in T. A simple calculation using Taylor's 
Theorem shows that 


2 3 

M 21 = ~ ^e* cos2 ( t n ^ + 0(/JT“ ) f° r a H (k,e) € T. 


Clearly, is not uniformly bounded in T for any choice of the constants 

* 

a >0 and kg as e 0. Consequently, the one-step method (5.6) is not 
uniformly (k ,e )-bounded in T. By Part b) of Lemma 3.1 we conclude that the 
one-step method (5.6) is not uniformly (k,e)-stable in r for any choice of 

k 

the constants a and kg. 

We now examine some numerical results which illustrate the practical 

T 

importance of uniform (k,e)-stability. Let v(t^) = (v (t ) » v n(t n ) ) and 
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v S (t ) = fv^(t ),v^ T (t denote the exact solution of (5.1) and the slowly 

varying part of v(t n ), respectively, at the grid points t n = nk, 

n = 0,1,***, with k = .1. Let x n = (x^,x^-|.) T and y n = (y™>yji) T denote 

the solutions of (5.5) and (5.6), respectively, at the grid points t n = nk. 

v(t n ), v s (t n ), x n an d y n were computed in double precision on the IBM 3081 

at Brown University. Values of these quantities at selected grid points are 

-3 

listed in the following tables for the three cases e = 2.5 x 10 , 

e = 2.5 x 10 -6 and e = 2.5 x 10 _9 . 


TABLE 1 
e = 2.5 x 10 


n 



n 

x i 

n 

y i 

0 

.100D+01 

.251D+02 

.100D+01 

. 100D+01 

1 

-.881D-01 

-.881D-01 

-.374D+00 

-.459D+00 

10 

-.308D+00 

-.308D+00 

-.294D+00 

-.525D+00 

50 

.639D-02 

.639D-02 

.608D-02 

. 109D-01 

75 

-.509D-03 

-. 509D-03 

-.486D-03 

-.868D-03 

100 

.240D-04 

.240D-04 

.230D-04 

.410D-04 

114 

. 100D-04 

. 100D-04 

.956D-05 

. 171D-04 

n 

v II (t n ) 


n 

X II 

a 

y n 

0 

. 100D+01 

. 100D+01 

. 100D+01 

. 100D+01 

1 

.900D+00 

.9000+00 

.813D+00 

.1710+01 

10 

. 199D+00 

.1990+00 

. 189D+00 

.339D+00 

50 

. 187D-02 

. 187D-02 

. 180D-02 

.319D-02 

75 

.189D-03 

. 189D-03 

.180D-03 

.323D-03 

100 

-.372D-04 

-.3720-04 

-.360D-04 

-.636D-04 

114 

.4260-05 

.426D-05 

.408D-05 

.727D-05 




















n 



n 

X I 

0 

.100D+01 

.250D-05 

. 100D+01 

1 

-.903D-01 

-.903D-01 

-.419D+00 

10 

-.310D+00 

-.310D+00 

-.295D+00 

50 

.646D-02 

.646D-02 

.616D-02 

75 

-.519D-03 

-.519D-03 

-.496D-03 

100 

.247D-04 

.247D-04 

. 237D-04 

114 

.103D-04 

. 103D-04 

.984D-05 

n 

V II (t n ) 

V n (t n> 

n 

X II 

0 

•lOOD+Ol 

. 100D+01 

.lOOD+Ol 

1 

.900D+00 

.900D+00 

.807D+00 

10 

.199D+00 

. 199D+00 

. 189D+00 

50 

. 191D-02 

. 191D-02 

. 184D-02 

75 

. 192D-03 

. 192D-03 

. 182D-03 

100 

381D-04 

-.381D-04 

-.364D-04 

114 

.441D-05 

•441D-05 

.4240-05 


. 100D+01 
. 931D+02 
.232D+03 
.486D+01 
. 391D+00 
. 186D-01 
.777D-02 


. 100D+01 
.925D+03 
.1496+03 
. 144D+01 
. 144D+00 
.287D-01 
.332D-02 























TABLE 3 


e = 2.5 x 10 


n 

V‘„> 

v®(t ) 
I n 

n 

x i 

n 

y I 

0 

. 100D+01 

.250D-08 

. 100D+01 

. 100D+01 

1 

-.903D-01 

-.903D-01 

-.419D+00 

-.927D+05 

10 

-.310D+00 

-.310D+00 

-.295D+00 

-.232D+06 

50 

.646D-02 

.646D-02 

.616D-02 

.486D+04 

75 

-.519D-03 

-.519D-03 

-.496D-03 

-.390D+03 

100 

.247D-04 

•247D-04 

.237D-04 

. 186D+02 

114 

. 103D-04 

. 103D-04 

.984D-05 

.776D+01 

n 

V II (t n ) 

V II (t n> 

n 

X II 

n 

y n 

0 

. 100D+01 

. 100D+01 

. 100D+01 

. 100D+01 

1 

.900D+00 

•900D+00 

.807D+00 

•924D+06 

10 

. 199D+00 

. 199D+00 

. 189D+00 

. 149D+06 

50 

. 191D-02 

. 191D-02 

. 184D-02 

. 144D+04 

75 

. 192D-03 

.1920-03 

. 182D-03 

. 144D+03 

100 

-.381D-04 

-.381D-04 

-.364D-04 

-.287D+02 

114 

.441D-05 

.4410-05 

.424D-05 

.3320+01 


Remark 5.1 . We list both the exact solution of (5.1) and the slowly 
varying part of the solution of (5.1) in order to illustrate Property 1.2. 

These numerical results permit us to draw the following conclusions for 
each value of e : 

1. After a transient phase, the solution of the one-leg method (5.5) 
approximates the slowly varying solution of (5.1) accurately to one 
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decimal place. The solution of the one-step method (5.6) does not 
accurately approximate the slowly varying solution of (5.1). 

2. The solution of the one-leg method remains bounded as e 0. The 

solution of the one-step method becomes unbounded as e + 0. 

According to calculations (5.8) and (5.9), the only way to reduce the 
size of the solution of the one-step method (and to accurately compute 
the solution or slowly varying solution of (5.1)) is to reduce the 
step size k. Consequently, the one-leg method is more efficient than 
the one-step method. 

Remark 5.2 . In order to efficiently solve stiff differential equations 
implicit methods are always used. The concept of uniform (k, e )-stability for 
system (5.1) or (1.4) permits us to compare the efficiency of implicit 

formulae. That is, an implicit method which is not a uniformly (k,e)-stable 
approximation of (5.1) ((1.4)) will be less efficient than an implicit method 
which is a uniformly (k,e)~stable approximation of (5.1) ((1.4)). 

The tables also illustrate conclusions 1-3 which are stated at the end 


of Section 3. 
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APPENDIX 


I. Notation 

I^ 3 * 5 * * * * * II denotes s-diraensional real Euclidean space, 
denotes s-dimensional complex space, 
denotes the transpose of the vector y. 
denotes the conjugate transpose of any matrix A. 
(i) 


<j: s 

T 


I! y II = max |y 
i 


denotes the vector norm. 


II A|| = sup II Ay II / II y [] denotes the subordinate matrix norm. 
y*0 

A matrix or vector f (t ,e) € C p (t,B) is all derivatives - — f ^ 1 ’ £ - with 

dt V 

0 < v < p are continuous and uniformly bounded by the constant B for 
0 < t < T and 0 e S. e g where > 0 Is a constant. 

A matrix or vector f (t,e) € C P (t,e,B) if all partial derivatives 


3 V f(t,e) 

v v 2 

at 3 e 


with 0 ^ + v 2 = V ^.P are continuous and uniformly bounded by the 

constant B for 0 t T and 0 _< e < Eq where Eq > 0 is a constant. 

Let f(E^, •••,£) and g(E^,»»*,e n ) be two vectors or matrices which 
depend on the n parameters anc * are defined for all 

contained in some set I. Then f(e^, •••,£) = 0(g(e^ , • • • , e^) ) for all 

(e i , • • • € I (or just f(e^,»**,e n ) = 0(g(ej , • • • ,e n ) ) when the set I is 
understood) if there exists a constant K, independent of £p •••»£> so that 

II f (e A • , e n ) II <_ K||g(e 1 ,..* ,£ n )ll for all (e^ , • • • ,e n ) € I. 
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II. Proof of Lemma 2.2 

We present the proof of Lemma 2.2 in this part of the appendix. I would 
like to thank Professor H.-O. Kreiss for helpful discussions pertaining to 
this proof. 

Proof of Lemma 2.2 . (All constants which appear in this proof are 
independent of k, t and e.) 

Write (2.8) as its equivalent one-step method 

(A.l) w n = a(k,t n> e)w n \ n = r, r+1 , • • • , and all (k,e)£ T' 

, n , n n-1 n-r+KT , 

where w = (y ,y , ..«,y ) and 


(A. 2) 


a(k,t,e) 


Q 1 (k,e,x 0 (t) 

I 

0 

• 

0 


T r (t))***«Q r (k,e,T 0 (t),...,T r (t)) 

0 • • • • 

io.. 

.... . 

...I 0 


Assumption (2.9) implies that a(k,t,e) is uniformly bounded for all 
rk t < «> and (k,e) £ T ' . 

Assumptions (2.9) and (2.10) imply that for all rk t < « and 
(k,e) £ r' 

(A. 3) o(k,t,e) = 8(k,t,e) + ky(k,t,e) 


with 
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(A.4) 8(k,t,e) = 


Q 1 (k,e,T 0 (t),...,T 0 (t))...*Q r (k,e,T 0 (t),.*.,T 0 (t)) 


I 

0 


0 • • • 

I... 


and 


(A. 5) 


y(k,t,e) 


E , • • • • E 

1 r 

0 o... 

0 0 • • • • . 


0 ...0 0 


Furthermore, g(k,t,e) and yCk,t,e) are uniformly bounded for all 
rk ^ t < » and (k,e) €T'. 

Let K^(k,t,e) and A^(k,t,e), i = l,»»*,mr, denote the eigenvalues of 
8(k,t,e) and a(k,t,e), respectively. By Assumption 2) (lines (2.11) and 
(2.12)) the eigenvalues of g(k,t,e) satisfy 

|K i (k,t,e)j k < 1 for i = l,««»,mr, 

(A. 6) and all rk _< t < °° and (k,e) € T' 

where k is a constant. The matrix y(k,t,e) is uniformly bounded, so there 

A A 

exists a constant k staisfying 0 < k <. F and a constant 0 _< X < 1 such 
that 

|x i (k,t,e)| < 1 for i = !,•••, mr. 


(A. 7) 


and all rk < t < 


00 


and (k,e) £ r" 
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where r" = ((k,e)£l? : 0 < e £q and a*e <_ k <. k} . 

& 

We now show that for each point rk < t < “ and any point (k,e) € T" 
there exist constants K and t, independent of t , k and e, and an 
invertible matrix S(k,t ,e) satisfying 

^ — 1 Vt 

(A. 8) sup { II S ( k , t ,e)ll, I! S (k,t , e ) II } <_ K < °° 

0_<t <» 

(k,e)er" 


such that \ 

(A. 9) II S *(k,t ,e)a(k,t,e)S(k,t , e ) II _< (x + ^ — ) < 1 for t e [t ,t + x] 

and all (k ,e ) € r" . 

By Shuras Theorem, for each rk <_ t < »o and (k,e) e r" there exists a 
unitary transformation U(k,t,e) such that 

(A. 10) U*(k,t,e)a(k,t,e)U(k,t,e) = a(k,t,e) = D(k,t,e) + P(k,t,e) 

where 

(A. 11) D(k,t,e) = diag(x i (k,t,e) , • • • ,X mr (k,t,e)) 

and P(k,t,e) is an upper triangular matrix with zeros on its diagonal. The 
matrices oi(k,t,e) and D(k,t,e) are uniformly bounded, so 

(A. 12) sup |(p(k,t,e))..|_<P<” for all 1 _< i, j _< mr 

Oj<t<oo ±J 

(k,e)er" 


where P is some positive constant 
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Let 

(A. 13) S - diag(l,5,...,5 mr “ 1 ) g = 4(mr-lTp ‘ 

Then straightforward estimates show that 

(A. 14) || (S) Lckjt.eJSn £ 1 ^ ^ for all rk jC t < “ and (k,e) € r". 

* 

By Assumption (2.18), for any rk t < « 

& & 

(A. 15) ct(k,t,e) = a(k,t ,e) + a(k,t,e) - a(k,t ,e) 

with 

(A. 16) |la(k,t,e) - a(k,t ,e)ll £ Q(t - t ) for all (k,e) € r". 

Set 

(A. 17) S(k,t*,e) = U(k,t*,e)l3. 

tfc 1 — inir 

Then S(k,t ,e) satisfies (A. 8) with K = £ . Furthermore, if we set 

x =— — , then straightforward estimates using (A. 15) - (A. 17) establish 

4QK 

(A. 9) and the preliminary result is justified. 

For our next preliminary result we set 

(A. 18) T* = min{x,k} and m = X + ^ ^ — < 1. 

Let q be the smallest integer satisfying u)^K < 1 and set 
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(A. 19) 6 = ojK 1/q and k A = — . 

0 q 

(Clearly 0 < 6 < 1.) We now show that for any integer £ > r-1 and any set 
of q+1 consecutive grid points t = (£+i)k, i = 0,»««,q, the solution of 
(A. 1 ) satisfies the estimate 


(A. 20) 


II w 


l+i„ 


< Ko) i ||w £ || 


for i = 0, 


, q and all (k,e) £ F 


where T is defined on line (1.13) with kg defined on line (A. 19). 
Consequently, by the first relation on line (A. 19) 


(A. 21) 


llw £+q i| < (6) q llw £ || for all (k,e)€T. 


To justify the claim, make the change of variables 


(A. 22) 


w £+i = S(k,t^,e)v £+i , i = 0 , • • • , q , 


in equation (A.l) where S(k,t,e) is the matrix defined on line (A. 17). 
(Note that S(k,t^,e) is constant over the (q+1) grid points.) Then 


(A.23 ) v £+i = (S 1 (k,t £ ,e)a(k,t £+j ,e)S(k,t £ ,e))v £+i 1 , i = l,...,q. 


By result (A. 9) and the definitions on line (A. 18) 


0+i i Q 

II v II <io P v II for i = !,•••, 


(A. 24) 


q, and all (k,e) £ T 
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By relation (A. 22) and (A. 8) 

(A. 25) ||w £+i || <_ KojSw^II , for i = l,»»*,q, and all (k,e) £ r 

and estimates (A. 20) and (A. 21) are justified. 

To complete the proof of Lemma 2.2, consider equation (A.l) and set 
n-j = £q+P where £ and p are non-negative integers with 
0 < p < q-1. By results (A. 20) and (A. 21) 

,w n . - 1 »*'> + p + h < K„ p .„ t ’ + h < 

for all (k,e)€r. Now set 6 = max{w,S} <1 to obtain 

II w n || < K<S n ^ H w^ II , for all r-1 < j < n < =° and (k,e) £ r. 


The proof of Lemma 2.2 is complete 
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